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S ! Abstract 

We revisit the proposed theoretical model for a small but nonzero cos- 
mological constant which seems supported increasingly better by recent 
observations. The model features two scalar fields which interact with 



each other through a specifically chosen nonlinear potential. We find a 
very sensitive dependence of the solutions of the scalar field equations 
on the initial values. We discuss how the behavior is similar to and dif- 
\Q • ferent from those in well-known chaotic systems, coming to suggest an 

interesting new type of the dissipative structure. 

The cosmological constant problem has two faces; an upper bound and a lower bound. 
q-i The former is given by the critical density p cr , or Q\ = A/p CT <l. In the Planckian unit 

system with 87rG = c = h — 1, we know p cr = 3Hq ~ t^ 2 , and the present age of the 
universe to ~ lOGy is of the order of 10 60 ; hence the upper bound A o b<,10~ 120 . On the 
other hand, almost any theoretical models of unification gives A t h ~ 1. The discrepancy 
is the well-known disaster. Recent observations seem to indicate a lower bound as well; 
Qa ~ 0.7. Let us try to understand this "small but nonzero cosmological constant." 

A solution for the upper bound may come from a scalar field <fi of the Brans-Dicke 
type: 

C = V=g ^<I> 2 R - \tg^d^d v <P + A + L mattcr ) , (1) 

where £ is a dimensionless constant related to BD's parameter u by = 1/4, and e = ±1. 
The scalar field contribution is shown to cancel the effect of A asymptotically. To see this, 
it is most convenient to apply a conformal transformation moving to the Einstein frame 
(as denoted by *) with nonminimal coupling removed. The same Lagrangian as Eq. |I] is 
now put into 

£ = V^gl UR* - \g^d, L ad u a - V(a) + L* mattcr ) , (2) 



*Based on the talk delivered at XXXIIIrd Rencontres de Moriond, Fundamental Parameters in Cos- 
mology (Les Arcs, France, January 17-24, 1998) 
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where the new canonical scalar field a is related to the original 0; <p = £ 1 ^ 2 e^ cr , with 
( = (6 + e^ 1 )- 1 / 2 . We also find that A now acts as a potential V(cr): 

7(a) = Ae" 4 ^, (3) 

which pushes a toward infinity; the effect of A continues to decrease. 
In spatially flat RW cosmology, the cosmological equations are 

3H 2 = p a + p m = p, and a + 3H& + V'{a) = 0, (4) 

where we omit the symbol * to simplify the notation. We should also supply the equation 
for the matter density p m as usual. Notice that p a = \b~ 2 + V(o~) is the energy density of 
the scalar field to be interpreted as an "effective cosmological constant" A e g. 
As an asymptotic solution, we obtain (for assumed radiation dominance) 

a(t) = t 1/2 , and a(t) = a + ^C 1 lnt, (5) 
Pit) = ^(l-^- 2 )*- 2 , with p„(t) = ^C 2 t- 2 . (6) 

The last equation shows that A c fj = p a does decay like t~ 2 , hence implementing the sce- 
nario of a "decaying cosmological constant." Today's A is small only because our universe 
is old 0. 

We notice, however, that A e g(t) falls off in the same way as the ordinary matter density 
Pmatter(^) does. Not qualified to be called a "constant," it is simply another form of dark 
matter. It does not help understanding the lower bound. We must include something that 
behaves more or less like a constant or at least falls off more slowly than t~ 2 . On the other 
hand, the behavior ~ t~ 2 was the key to have a "small" A. As a compromise we expect 
that A e g(t) would decay like ~ t~ 2 as an overall behavior, but with some "landings," in 
one of which we today happen to be, and which mimics a cosmological "constant" for 
some duration of time. 

Trying to implement the idea, we came across a model on a try-and-error-basis, by 
introducing another scalar field $ which couples to a |||3|]. The Lagrangian is the same 
as Eq. |] with the added kinetic term of $ and V(o~) replaced by (see Fig. 2 of Ref. 2) 

V(a, $) = e- 4CCT + ^m 2 $ 2 [1 + p sin(K(r)A , (7) 

where m, [3, k are the constants naturally of the order unity. It has a "central valley" 
given by Eq. ^, but shows an oscillatory behavior ~ sin(/ter) if we climb the wall in the 
direction of 
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A typical solution is shown in Fig. |1|. We started the integration at some time in the 
post-inflation era. We also mildly fine-tuned initial conditions such that somewhat 
smaller than 1, say, at to ~ 10 60 . Notice a sequence of alternate occurrence of a rapid 
change and a nearly standstill of the scalar fields. Also toward the end of each "landing" of 
p s (= A e g(t)), which is now the energy density of the coupled cr-$ system, the scale factor 
a(t) shows a mini-inflation, an expected behavior. Magnifying a portion of Fig. |l](b) 
around the present time (log to ~ 60) extending to a tiny interval (Alogt ~ 1) yields a 
"practical" plot going back to z ~ 5, showing A e g- which looks nearly constant. 

These landings, as indicated also by other examples, tend to occur nearly periodically 
with respect to Int. This seems to suggest that some cosmological anomalies might have 
taken place nearly periodically with respect to Int. 

As we find, the results depend on the choice of the parameters quite sensitively, par- 
ticularly on the initial conditions. In this connection, we first point out that the most 
important ingredients are present in the dynamics of the system of the two scalar fields. 
We find the alternate occurrence of the "wake-up phase" and the "dormant phase" of the 
scalar fields even without the cosmological environments, provided the frictional coeffi- 
cients 3H vary as ~ t . 

Fig. ^ is one of the examples of the solutions of the isolated 2-scalar system with the 
same potential with the frictional constants T = j/t. Suppose we start from an initial 
position off the central valley on the potential Eq. [7[ o will be pushed in the positive 
direction. The potential will decrease quickly because it falls off like e~^ a . Soon both 
fields will be virtually free. Due to the frictional forces they will be decelerated almost 
to a complete stop. However, as the time goes on, the frictional forces will also decrease 
according to ~ t _1 . They will become even weaker than the forces derived from the 
potential which had much dwindled. Then the fields are suddenly pushed again. But the 
potential will be weakened soon again leading to the next dormant period. 

However, in order for the fields to be pushed strongly after the dormant phase, their 
positions and velocities must "match" the previous values to some extent. Otherwise, the 
fields will not be waken up, beginning simply to fall the potential slope slowly toward 
and along the central valley. This is the asymptotic behavior, as we find in this example, 
reached commonly for any initial conditions. Obviously this corresponds to a "fixed-point 
attractor," like in standard damped oscillators. 

The "matching condition" is rather subtle. This is the reason why a slight change 
of initial values may sometimes result in a considerable change in the behaviors after- 
ward. This sensitivity to the initial values reminds us of a well-known chaotic behavior 
in nonlinear systems. Sudden change of the number of repetitions before reaching the 



3 



smooth asymptotic behavior, for example, might be compared to sudden occurrence of 
"bifurcation" in typical chaotic solutions. 

However, there is an important difference from purely chaotic systems. Usually a 
chaotic behavior is characterized by "strange attractors," while our solutions tend to a 
simple fixed-point attractor. Complicated hard-to-predict states occur only during a finite 
transient period. In accordance with this no fractal structure is produced in the orbits in 
phase space; fractal structure might be present only to a limited "depth." 

We thus find what might be called an "incomplete chaos" or an "incomplete fractals." 
They can be still interesting because genuine chaos might be a mathematical idealization 
for natural phenomena. Also chaotic behaviors are seen almost every aspect of life. But 
obviously life lasts only during a limited lifetime; the ultimate destination is a non-chaotic 
death, corresponding to a fixed-point attractor. In this sense, our model might provide 
a primitive example of a living system. On the other hand, it seems to give a model of 
pattern formation (with variable period). 

In our model, near periodicity in In t is a consequence of T ~ t~ x . There are many 
examples of pattern formation with some periodicity. Patterns on sea shells, for example, 
can be reproduced by certain nonlinear equations of rather simple structure [|]]. Periods 
of these patterns are not determined by the restoring force but, unlike in simple systems, 
primarily by friction or dissipation in addition to the coupling strength. This feature is 
shared in our model. In this sense also, our nonlinear equations might provide another ex- 
ample of "dissipative structure." It is interesting to notice that a model of game dynamics 
with more than three species also generates periodicity with respect to hit |J. 

Did we introduce as many degrees of freedom and parameters as we need? The answer: 
the cosmological constant problem is so challenging that we may not have reached even 
close to any promising result no matter how many ingredients we introduced if we have 
come in a wrong direction. It seems we are on a right track. 

We thank Kunihiko Kaneko and Takashi Ikegami for enlightening discussions on non- 
linear systems. 
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Figure 1: (a) In a (solid; a for the scale factor), a (dotted), 2$ (broken) against logi. 
(b) logp s (solid), logp m (dotted). The parameters are A = 1,( = 1-582, m = 4.75, (3 = 
0.8, k = 10. The initial values at t x = 10 10 are a l = 6.75442, b x = 0, $i = 0.212, $ x = 0. 
At log t = 60.15 corresponding to the present age chosen to be 1.21 x 10 10 y, we obtain 
H = 81km/s/Mpc and f2 A = 0.67. See also Captions to Figs. 1 and 2 in Ref. 3. 
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Figure 2: An example of the two-scalar system, in which the solution enters the asymptotic 
behavior after two repeated patterns, a (solid) and 2$ (dotted) and the broken line for 
the asymptotic behavior of a. $ will tend to zero slowly. The parameters and the initial 
values are the same as those in Fig. 1 except for the frictional coefficient chosen to be 
(3/2)t _1 . See also Caption to Fig. 3 in Ref. 3. A slight change of the initial values may 
result in a sudden change of the repetition number to three, for example, also changing 
the time of the end of the repetition period. 
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